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SUMMARY 


The development and the numerical application are presented of a 
Rayleigh-Ritz, or modal, type of flutter analysis which takes into 
account three-dimensional structural and aerodynamic behavior. The 
flutter mode is approximated by a series of natural -vibration modes, 
and the aerodynamic forces corresponding to these modes are derived 
from subsonic lifting-surface theory, according to the kernel-function 
approach, for a finite wing oscillating in compressible flow. 

The application is made to a delta semispan wing with a leading- 
edge sweep angle of 45° which fluttered at a Mach number of O.85. Results 
of flutter calculations show that, for this case, when the first three or 
four natural -vibration modes are used to approximate the flutter mode, 
converged solutions for the flutter speed are obtained that are about 
5 percent less than the experimental value. Theoretical flutter-speed 
boundaries were located for a range of densities and Mach numbers 
including those of the experimental-flutter condition. Further applica- 
tion of the analysis to study the effects of variation in certain struc- 
tural properties showed that the converged flutter speeds were more sen- 
sitive to variations in the natural frequencies than to either variations 
in mass or to the inclusion of generalized-mass coupling terms whose 
existence is due to the use of experimental natural mode shapes. 


INTRODUCTION 


Current aircraft -design trends, such as the use of thinner struc- 
tures and external stores on aircraft capable of very high speeds, have 
combined to diminish flutter safety margins and have, consequently, 
increased the need for greater acciiracy in flutter prediction. As a 
result, both the structural and the aerodynamic aspects of the flutter 
problem should be treated by use of more realistic methods than by the 
beam- theory and strip-theory methods commonly employed in the past. 
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This report illustrates the development and application of a method 
of flutter analysis which takes into account three-dimensional structural 
and aerodynamic behavior. In treating the structural problem, the flutter 
mode may be approximated by a series of either natural or assumed vibra- 
tion modes which could have platelike distortions or shapes. In treating 
the aerodynamic problem, lifting-surface theory is used to obtain aero- 
dynamic forces which take into account finite span and compressibility 
as well as the modes of vibration of the structure. 

In order to illustrate the application of the method, a number of 
flutter calculations are performed. Primary attention is directed 
toward correlation of the calculated result with an experimental flutter 
result for a delta semispan wing with a leading-edge sweep angle of 45 ° 
at a Mach number of O.85. The calculations are based on carefully meas- 
ured natural -vibration modes obtained by means of an optical method. A 
detailed description of the determination of these modes and of the mass 
distribution is given in appendix A. 

Other calculated results are presented for the same wing plan form 
in order to show some effects of variations in air density and in Mach 
number (for a Mach number range from 0 to O.95). Related questions con- 
cerning the mmiber of modes required for convergence and certain effects 
of nonorthogonality of the measured modes are considered. N\jmerical 
evaluation of the elements of the flutter determinant is discussed in 
appendix B. 


SYMBOIS 


( j) 

b 

d 

Fn(6) 

f 


generalized aerodynamic force (see eq. (8)) 

arbitrary constant in series form of pressure dis- 
tribution (see eq, (l?)) 

local wing semichord, ft 

wing root semichord, ft 

distance between mirrors on wing surface and screen 
for optical method of measuring modes, in. 

fimction based on chordwise pressure term (see 
eq. (B 2 )) 

circular frequency, cps 
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^’nm 

g 

h(x,y,t) 

hi 

Ic 

Is 

^nm 


K(M,ko,x-r,y-ii) 

Lj 

I 

M 

Mi 


Mij 

m(x,y) 

N 


constants used in numerical integration (see 
table VI) 

dimensionless pressure mode (see eqs. ( 17 ) ( 18 )) 

structural damping coefficient 

instantaneous deflection of point on wing surface 
in the flutter mode , ft 

normalized displacement at point x,y in ith mode of 
vibration^ hi(x,y) 

chordwise integrating factor appropriate to station c 
(see eq. (B5)) 

spanwise integrating factor appropriate to station s 
(see eq. (b6)) 


surface integral in generalized aerodynamic force 
(see eq. (B2)) 

kernel function of three-dimensional integral equa- 
tion (see eq. (l2)), l/sq ft 

dimensionless form of three-dimensional kernel func- 
tion (see eq. (l4)) 

reduced-frequency parameter , bgCD/V 

dimensionless pressure function (see eq. ( 5 )) 

wing semi span, ft 

Mach number 

generalized mafes associated with ith mode of vibra- 
tion (see eq. (3))^ slugs 

generalized-mass coupling term (see eq. (2l)), slugs 

local mass per unit area at point x,y, slugs/sq ft 

nximber of sheets of aliiminum foil on top or bottom 
surface of wing 
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Ap(x,y,t) 


Apj 


q 


i 


S 

t 


local oscillating pressure difference between top 
and bottom surfaces of wing in the flutter 
mode (see eq. (^)), Ib/sq ft 

pressure difference at point x^y in jth mode of 
oscillation^ Apj(x^y ) 

generalized coordinate in ith mode of oscillation, 


complex amplitude of generalized coordinate in 
ith mode , ft 

area of wing surface, sq ft 
time, sec 


ti,tjp,t£ 


thicknesses of aluminTom insert, aluminixm foil, and 
balsa, respectively (see eq. (AD) 


V 




Z 

a 

a, a' 




velocity of airstream, fps 

weight of wing per unit area, Ib/sq in. 

Cartesian coordinates (see sketch following eq. (l2)) 

dimensionless chordwise variables referred to bQ 
(see eq. (l4)) 

dimensionless spanwise variables referred to Z 
(see eq. (ih)) 

airfoil ordinate (see eqs. (A2) and (A5)) 

local slope of wing in pitch during vibration, Sh/Sx 

local slope of wing in a direction normal to line 
(or ray) of constant percent chord 

vectors representing angular displacements a and 
, respectively, according to the right-hand 
vector rule (see eqs. (A 5 ) and (a 6 ) and fig. 7 ) 

unit weights of aluminum alloy and laminated balsa, 
respectively (see eq. (A1)), Ib/cu in. 

reflected displacement on screen in optical method 
for measuring mode shapes, in. 
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5 


6x,6y 


streamwise and spanwise components, respectively, 
of 6 


0 


A 


r 


l^i 


0 

0 ’ 


0 , 0 ' 

n 


cu 


Indices : 
c, s 
f 


angular chordwise variable (see eq. (19)), deg 

angle of sweepback for rth ray, deg 

generalized mass-density ratio in ith mode of vibra- 
tion (see eq. (lO)) 

air density, slugs/cu ft 

local slope of wing in roll during vibration, Sh/Sy 

local slope of wing in a direction parallel to line 
(or ray) of constant percent chord 

vectors representing angular displacements and 0 ' , 

respectively, according to the right-hand vector rule 

angle of image, referred to the horizontal, on screen 
in optical method for measuring mode shapes, deg 

complex eigenvalue of determinantal flutter equation, 

^ is) 

angular frequency^ 2 itf^ radians /sec 


chordwise and spanwise stations^ respectively 
flutter 


J 


modes of vibration under consideration 


n,m 


r 


X 


y 


chordwise and spanwise pressure modes^ respectively^ 
in aerodynamic quantities (m denotes power to 
which spanwise variable is raisedj see eq. (l 8 )) 

line (or ray) of constant percent chord on wing 

streamwise direction (related to local pitching 
direction) 

spanwise direction (related to local rolling 
direction) 
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Matrix notations: 



rectangular matrix 
row matrix 

column matrix 
diagonal matrix 

Dots over symbols denote differentiation with respect to time. 


METHOD OF FLUTTER ANALYSIS 


The present section is concerned with the development of a working 
form of the flutter equations of a Rayleigh-Ritz, or modal, type of 
analysis based on aerodynamic forces obtained by subsonic lifting-surface 
theory. A brief discussion of the method of obtaining the aerodynamic 
forces is included. 


Development of Flutter Equations 

In the development of the flutter equations, a basic assuir^jtion of 
the Rayleigh-Ritz approach is that the displacement h(x,y,t) corre- 
sponding to the flutter mode may be represented by a superposition of 
either natural or assumed modal fimctions in the form 


h(x,y,t) = q2^hj_ + q2h2 + . . . (l) 

where qj^ = qj^e^^ is the generalized coordinate of the ith chosen mode 
and hj[ = hi(x,y) is the corresponding normalized mode shape. 

In the present Investigation the normalized modal functions hi 
are chosen as the natural (orthogonal) modes of vibration of the struc- 
tiire under consideration. With this choice of mode shapes, the general 
equation of motion (obtained, for exan^jle, as in ref. l) in the ith 
degree of freedom may be written in the form 


Mill + ^^^Mi^i 



x,y,t) 


dS 


• ( 2 ) 
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where Mi represents the generalized mass in the ith mode and is defined 
as 


Mi = 



(5) 


The integral on the right side of equation (2) represents the generalized 
aerodynamic force and contains the aerodynamic loading Ap(x,y,t) which, 
consistent with equation (l), is regarded as a superposition of 
aerodynamic -loading modes 


-^pCxjyjt) = q^_ Apj^(x,y) + qg Ap2(x,y) + . . . (4) 


where Api:(x,y) denotes the aerodynamic loading associated with the 
J 

mode shape hj , In application^ with 
expressed in terms of dimensionless functions 


— io3t 




^(x,y,t) is 
through the relation 


Ap(x,y,t) 



qp 

+ — Lp + 
^0 



(5) 


The form of the complex functions Lj is dealt with at a later stage. 
At this stage, use may be made of equations (l) and (5) to obtain equa- 
tion (2) in the form 


1 



4jtpV^Z 

bo^o]^ 


jj' + 

S 


qpLp + 


.)dS = 0 (6) 


or, alternatively. 



1 


+ • 


. .) = 0 


(7) 
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where 


If ^ 

S 

and 

= ^ 


If equation (7) is written for each of the degrees of freedom under 
consideration, a set of simultaneous homogeneous equations results. The 
flutter condition is then given hy the vanishing of the determinant of 
the coefficients of q^ in these equations; thus. 


1 



^11 


, 2 
^0 M-i 


^21 


1 


^12 

ko^^Al 


A22 

— a + 


. . . = 0 ( 9 ) 



where 


2^ 


4jtpl^bQ 

"mT 


(10) 


and where a is defined hy 


a = 



ig) 


( 11 ) 


In this form ft is a complex eigenvalue and contains the unknown fre- 
quency CD and a damping coefficient g that becomes zero at the border- 
line flutter condition between damped and undamped motion. 
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Equation (9), together with the definitions given in equations ( 3 ) > 
(8), (10), and (11); constitutes the vorking form of the flutter equa- 
tion. In other sections of the report the determination of the ingre- 
dients of equation (9) is discussed; appendix A is concerned with the 
determination of mode shapes and mass distribution for a particular con- 
figuration, and appendix B deals with numerical techniques used to 
evaluate generalized masses and generalized aerodynamic forces. 


Determination of Aerodynamic Loading 

In the determination of the aerodynamic loading, the functions Lj 

required in the generalized aerodynamic forces are obtained from the 
integral equation which relates lift and downwash distributions in sub- 
sonic lifting-siirface theory. (See, for example, ref. 2 .) A systematic 
numerical solution is employed herein which is similar to that of refer- 
ence 5 but which makes use of a more exact form of the kernel function 
together with more refined numerical-integration procedures. The method 
employed has been programed for the IBM type 704 electronic data pro- 
cessing machine . 

For the purpose of the present investigation, the integral equa- 
tion may be written as 

-(s " t) 

^ s 

where the coordinate-axis system is shown in the following sketch: 
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where ^te <ienote, respectively, the coordinates of the 

leading and trailing edges. The term on the left side of equation (l2) 
denotes the downvash angle at point x,y associated with a displacement 
h(x,y,t)* On the right side of equation (l2), the kernel function 

K^M,^,x-^,y-T)j represents the downwash produced at point x,y by a unit 

pressure load at point the function Ap(|,T],t) represents the load 

required to satisfy the specified downwash condition. (The minus sign 
on the left side of eq. (l2) is associated with the form of the kernel 
function and arises from the use herein of the sign convention for down- 
wash, displacement, and lift which are positive downward.) 

Substituting equations (l) and (5) into equation (l2) gives 



Introducing the following dimensionless quantities 


X 










> 




( 1 ^) 
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into equation (l3) gives 


A + ikoV ^ hi + ^ h2 + 


idx 


'0 


'0 


K(M,ko,x-|,y-Ti) Li + — L2 + • 
’ ) ^0 


d| dq 


(15) 


From equation (l5) the form of the integral equation for the jth 
mode of oscillation may be seen to be 


+ ikojhj = j j K(M,ko,x-|,y-T]) Lj d? dii 


(16) 


In the solution of equation (l6), it is assumed that a function Lj may 
be represented by a series of pressure modes of the form 


n m 


^nm 


(17) 


where the quantities a^^ are arbitrary constants to be determined. 

The indices n and m are associated with chordwlse and spanwise pres- 
sure modes, respectively. The fonu of the functions fj^j^ is dictated 

by known leading- and trailing-edge conditions. In subsonic flow they 
should satisfy the edge conditions pertinent to the Kutta condition. In 
the existing procedure for solving the Integral equation they have been 
expressed (in terms of an angular chordwise variable 0) as 


^Qm 


^'im 

^2m 


m 


-^br 


q \|1 - TJ cot I 
b ^ 


_m 


T 


sin 0 




sin 20 


(18) 
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where j for example, for symmetric motion the index m is equal to 
0, 2, 4, . . . and where the variable 9 (see sketch which follows 
eq. ^2)) is related to | by the expression 


I 


^Ze 


— cos 0 
^0 


(19) 


in which denotes the equation of the wing midchord. 

2 

If the series form of Lj given by equation (17) is substituted 

into equation (l6)^ a working form of the integral equation is obtained 
as 





(j) 

Q-nm 


n m 


fl 


^'nm K(M,kQ,x-^,y-il) d| dq 


( 20 ) 


and may be seen to consist of a summation of definite surface integrals^ 
each weighted by an iinknown constant - 


In order to determine the constants Q- collocation procedure 

is used. The right side of equation (20) is evaluated for as many 
points x^y^ designated as control points^ as there are unknown 
constants and is equated to the known downwash angle at each con- 

trol point. A set of simultaneous equations is thus obtained which may 
be solved for the values of armi- Once determined, the constants a^m 
are used with equations (17) and (18) to define the pressure function Ljj 
and Lj may then be employed in equation (8) to obtain the generalized 
aerodynamic forces Aj_j . In the present application nine terms of the 

series in equation (17) were used (so that n = 0, 1, 2 and m = 0, 2, 4), 
and the downwash was satisfied at the nine control points shown in fig- 
ure l(a). As previously noted, a solution of the integral equation as 
represented by equation (20) has recently been programed for the lEM type 
704 electronic data processing machine. 


APPLICATION OF METHOD TO A DELTA WING 


In order to show the application of the analytical treatment described 
in the preceding section, a number of calculations have been made for a 
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delta semispan wing with a leading -edge sweep angle of 45°. Primary 
attention is directed toward a correlation of calculated results with 
an experimental-flutter result. Other calculations are performed for 
the configuration of the experiment to explore some effects of varying 
Mach number and air density. In addition, some effects of varying mass 
and natural frequencies and certain effects of nonorthogonality of the 
measured modes are examined. 


Correlation of Calculations With Experiment 

The experimental result used as a basis for the calculations has 
been obtained by William T. Lauten, Jr., and Marvin F. Birrgess at the 
Langley Aeronautical Laboratory. Construction details of the 45° delta 
semi span wing under consideration are shown in figure 2. The mass prop- 
erties of the model are shown as chordwlse and spanwise distributions 
of weight per unit area in figure 3 . The mode shapes in the first four 
natural modes of vibration were obtained by means of an optical method 
and are shown in figure 4. Details of the methods used in obtaining 
mass and mode shapes are described in appendix A. 

Two models were involved in the experimental program. One, des- 
ignated as model A, fluttered to destruction before its mode shapes and 
mass distribution had been deterrtiined. These structural properties 
were obtained by use of a second model, model B, which was built to the 
same specifications and had very similar nodal patterns. As can be seen 
in table I, model B was lighter than model A and had different natural 
frequencies. In all flutter calculations the mode shapes for model B 
were used, and in the calculations for model A the generalized masses 
were adjusted by the ratio of the total masses of the models. 

The flutter calculation of primary interest has been made for 
model A, which fluttered as a cantilever at a Mach mmiber of O .85 with 
an air density of O.OOO 787 slugs/cu ft. Results of this flutter cal- 
culation are shown in table Il(a). Converged flutter-speed solutions 
are obtained (as indicated by the agreement of the first three- or fo\xr- 
mode calculations) which are about 5 percent less than the experimental 
value. Calculated results based on two two-mode subcases are included 
in table Il(a) . The flutter speed calculated with modes 1 and 5 is 
about 2 percent above the experimental flutter speed. It is noted that 
modes 1 and 5 resemble, respectively, first-bending and first-torsion 
modes of a beam. 


Variations in Mach Number and Air Density 

Results of converged flutter calculations at other values of Mach 
number and air density for the configuration of the experimental flutter 
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condition are given in table II (b) and are presented, in different forms, 
in figures 5 and 6. Figure 5 illustrates the effect on of varying 

density, with Mach number held constant at 0.85- Results are shown in 
terms of an air-density parameter l/^P* The effect of Mach mxmber on 
calculated flutter speed is shown in figure 6 for each of three values 
of air density p. Results are presented in terms of a flutter-speed 
coefficient Vf/bonj. As can be seen, the flutter speeds calculated for 

the present configuration appear to remain nearly constant for all Mach 
numbers up to M = 0.95 • The corresponding calculated flutter frequencies, 
on the other hand, tend to decrease with increasing Mach niaaber, as may 
be seen in table Il(b). 


Effects of Certain Structural Modifications 

The effects of the differences in total mass and natxiral frequencies 
between models A and B, noted in table I, were considered worthy of fur- 
ther study. First, flutter calculations were performed for model B, and 
the converged flutter speed was 10 percent less than that for model A. 
(Compare case 1 of table Il(c) with the four-mode result in table 11(a).) 
In an effort to separate the effect of the difference in total mass from 
the effect of the difference in natural frequencies, additional flutter 
calculations based on the generalized masses of model B and the natural 
frequencies of model A were performed. The results, listed in table Il(c) 
(case 2 ), show that about three-fourths of this 10-percent difference in 
flutter speed was due to the differences in the natural frequencies 
between the two models . The remaining difference in flutter speed is 
then attributed to the difference in total mass. 

As can be seen in equation ( 9 ), the offdiagonal elements of the 
flutter determinant contain no inertial-coupling terms. This is a con- 
sequence of the orthogonality condition for natural coupled modes, which 
may be expressed as 


M- 


ij 



h^hj dS = 0 


(1 / j) ( 21 ) 


It is recalled that, in the present analysis, measured modes were 
employed. In order to investigate their orthogonality, values of M^j 
were computed for model A by use of equation (21) and are listed in 
table III. In order to determine the effects of the nonzero values 
of they were added to the appropriate elements of equation ( 9 ) 

and two calculations were made for model A. First, the aerodynamic 
terms were eliminated and the natural frequencies were computed; results 
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are shown in table III and, with the exception of are seen to differ 

by a negligible amount from those in table I . In a second calculation the 
change in flutter speed due to including the terms was determined; 

the other parameters used were those of the four -mode solution of 
table 11 (a). Results are listed as case 3 of table Il(c) and show that 
the converged flutter speed has been decreased by less than 1 percent. 


CONCLUDING REIiARKS 


The development and the numerical application have been presented 
of a method of flutter analysis which takes into account three-dimensional 
structural and aerodynamic behavior. The flutter mode was approximated 
by a series of natural -vibration modes, and the aerodynamic forces cor- 
responding to these modes were derived from subsonic lifting-surface 
theory, according to the kernel-function approach, for a finite wing 
oscillating in compressible flow. 

The application was made to a delta semispan wing with a leading- 
edge sweep angle of 45° which fluttered at a Mach number of O. 85 . 

Results of flutter calculations show that when the first three or four 
natural-vibration modes were included, converged solutions for the flutter 
speed were obtained which were about 5 percent less than the experimental 
value. Further application of the analysis was made to orient the experi- 
mental flutter condition with theoretical flutter-speed boundaries for a 
range of densities and Mach numbers including those of the experiment and 
to study the effects of certain variations in natural frequencies and 
total mass . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., July 17, 1958- 
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APPENDIX A 

STRUCTURAL PROPERTIES OF DELTA-WING MODEL 


The structural properties of the delta-wing model needed for the 
flutter analysis that is presented and applied in this report consist 
of the mass distribution and the natural modes of vibration, which are 
represented in equations ( 5 ) and (8) as m(x,y) and h^, respectively. 

This appendix provides a description of the methods used to determine 
these properties for the delta-wing configuration considered in this 
analysis . 

Description of Model 

Figure 2 shows pertinent details of the construction of the model. 
This construction consisted of a 0.02-inch- thick sheet of aluminum alloy 
located at the plane of symmetry with vertically laminated balsa glued 
to each surface of the sheet. The balsa laminations were shaped to the 
ordinates of an NACA 65A004 airfoil section, and the outside surface was 
wrapped with layers of 0.001-inch-thick alvuninum foil. The wing tip was 
rounded so that the semispan was 52.85 Inches at a point I.5 inches for- 
ward of the trailing edge. The aerodynamic aspect ratio, with the assump- 
tion of a straight or squared-off tip based on extensions of leading and 
trailing edges to a common point, was 3-54. 

Two models were involved, both built to the same specifications. 

One model, designated as model A, fluttered to destruction before all of 
its properties were determined; the other model, model B, was used to 
obtain the mass distribution and experimental natural-mode shapes for 
these elements. Certain mechanical properties of both models are com- 
pared in table I . 


Mass Properties of Model 

The mass properties of both models studied in this Investigation 
are based on the chordwise and spanwise distributions of weight per unit 
area shown in figure 5* These distributions were calculated from the 
equation 

“v, - ’'a4i * 't) * Vb 

where w^ is the weight of the wing per vinit area, is the unit 

weight of alumimmi alloy, tj is the thickness of the aluminum insert. 



MCA TN 4595 


17 




tp is the thickness of the aluminum foil, tg is the thickness of the 
laminated balsa, and y-Q is the unit weight of laminated balsa. The 
value of which was determined from the measured total weight of the 

wing and the vol\mie occupied by the balsa, was found to be 0.01007 Ib/cu in. 
By using this value together with y^ = 0.100 Ib/cu in. and the data given 

in figure 2, equation (Al) can be written in the form 



everywhere on the wing except in the regions of the leading and trailing 
edges which were covered by an additional layer of aluminum foil which 

was li inches wide on each surface. The consequent increase in weight 

4 

in these edge regions was accounted for in the equation 

w = 0.100(0.022 + 0.002N) + 0.01007[48zb - (0.022 + 0.002N)] (A5) 

w ^ 

In equations (A2) and (A5) N is the number of 0.001-inch-thick 
sheets of aluminum foil on each (top or bottom) surface of the wing, and 
z is the ratio of one-half the local thickness to the local chord for 

an MCA 65A004 airfoil. 


Optical Method of Measuring Natural -Vibration Modes 

The natural -vibration modes were obtained from the results of shake 
tests conducted at zero airspeed and involving the use of an electro- 
dynamic shaker. The shaker was located near the wing root in order to 
reduce as much as possible the effect of the shaker mass on the natural 
modes. For each natural frequency found, the corresponding mode shape 
up through the fourth mode was determined by application of an optical 
method described in this section. This method involved the measurements 
and numerical integration of local slopes in pitch and roll directions 
at 24 stations distributed over the wing surface. 

A schematic diagram of the apparatus used in this optical method is 
shown in figure 7. Small mirrors that were l/l6 inch square were glued 
to the wing at six stations along rays at 25 , 50 , and 75 percent chord 
and along the trailing edge as shown in the mirror lattice in the upper 
left corner of the figure . Light aimed at the wing was reflected by the 
jj^pj*]^ors onto a screen located a considerable distance from the model. 
During vibration in a natural mode, images reflected from these mirrors 
appeared as straight lines whose lengths and directions were marked on 
the screen. These measured lengths and directions were resolved into 
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APPENDIX B 

EVALUATION OF ELEMENTS OF DETERMIMNTAL FLUTTER EQUATION 
AND APPLICATION TO A SPECIFIC CONFIGURATION 


Evaluation of Generalized Aerodynamic Force 


In order to evaluate the generalized aerodynamic force 

(eq. (8)), use is made of the dimensionless variables introduced by equa- 
tion (l4) and the pressure function L. given by equations (17) to (19) 

J 

to write equation (8) as 



In the evaluation of this equation, the following definitions are useful: 


F^(0) = cot — sin 0=1+ cos 0 
0 ^ ' 2 


Fj_(0) = sin 0 

F2(0) = ^ sin 20 sin 0 

1 pi 


= f f” iii 10 


0 JO 


(B2) 


-(i) 


Equation (B1) may be written in terms of the integral as the 

summation 





(B3) 


n m 
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or^ as the matrix product 



Numerical integrating techniques can he applied to evaluate the 
surface integrals in the row matrix of equation (b4) . (it is 


noted in this connection that the use of the angular coordinate 6 in 
the chordwise integration is particularly convenient for this applica- 
tion. For example^ when represented in Cartesian coordinates, the first 

term in the integrand of infinite at the wing leading edge; 

when expressed in polar coordinates, the product cot ^ sin 0 is a non- 
singular function. ) A number of spanwise stations may be chosen at 
y = with appropriate spanwise integrating factors Similarly, 

at each spanwise station, a number of chordwise stations are selected at 
0=0^ with appropriate chordwise integrating factors s and c 


thus* identify stations on an integrating lattice such as that shown in 
figure l(b). 
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The normalized displacements hp at the integrating stations may be 
arranged in a rectangular matrix of the form 



Equations (B 5 ) to (BT) niay be used to obtain the surface integrals 
of equation (B 2 ) from the matrix product 

m -4 

The elements of the rectangular matrix on the left side of equation (b 8 ) 

may then be rearranged as the row matrix of equation (b 4 ) and be 

used with the values of the constants a^*^^ to obtain the generalized 

, nm 

aerodynamic force A,- i . Equation (b 4 ) may, of course,, be expanded for 
any number of modes into 
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Evaluation of (Jeneralized Mass 

The generalized mass is expressed as a surface integral by 

equation (3) and may be evaluated by use of the same matrices of inte- 
grating factors and displacements employed in evaluating the generalized 
aerodynamic forces A. For this purpose^ equation (5) may be written 

in terns of the variables defined by equations (l 4 ) and (19) as 


i(01,y) ^ sin 0 d 0 dy 


(BIO) 


or^ as a summation which involves the integrating factors and 

used with A^^, in the form 

“l " “o* ^ X ^ ■"(<*c.ys)l'i^(ee,ys)=in 6^ (BH) 

s c 


or in terms of a matrix product as 



(B 12 ) 


The column matrix 



which pertains to a particular spanwise sta- 


tion^ may be obtained by squaring the elements in the corresponding 
column in the matrix ^hj (eq. (B 7 ))- The inertial-coupling term M^j 

(eq. (21)) is found by replacing the matrix in equa- 
tion (B 12 ) with a matrix of crossproducts for i ^ J . 


Application of Method to a Specific Configuration 

In the application to a specific configuration of the integrating 
procedures discussed in the previous section, either equal or unequal 
intervals in a given variable may be taken. With equal intervals, 
appropriate integrating factors are readily available; when unequal 
intervals are employed, corresponding integrating factors must be devel- 
oped. Useful procedures relating to numerical integration are presented 
in reference 5* 
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In the present application to the flutter of a delta wing, 10 equal 
intervals in the chordwise variable 6 and 8 equal intervals in the 
spanwise variable y were used. The integration stations are shown in 
figure 1(b) . Chordwise integrating factors employed are based on a 

rule employing overlapping quintic functions derived in reference 4. 
Spanwise integrating factors were based on a rule which employs 

overlapping quartic functions . Data pertinent to the evaluation of the 
factors and I„ are given in table VI. 

Normalized mode shapes employed in obtaining both generalized 

force Aj^j and generalized mass were obtained from the plots of 

figure 4. The normalizing station for each mode was arbitrarily located 
on the midchord at y = 0.875- Values of mass per unit area m(9c,yg) 
for use in determining were obtained from the weight distribution 

shown in figure 3- The components and hj of the downwash angle 

employed as a boundary condition in equation (20) are listed in table VII 
for the control points shown in figure l(a). 
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TABLE I 

CERTAIN MECHANICAL PROPERTIES OF DELTA-WING MODELS 


Item 

Model A 

Model B 

Aspect ratio (squared-off tip) 

Airfoil section 

5.5^ 

NACA 65AOO4 

5.5^ 

MCA 65 A 004 

deg: 

Leading edge 

25 percent chord (r = I) 

50 percent chord (r= 2 ) 

75 percent chord (r=3) 

Trailing edge (r= 4 ) 

^0" 

Z , f t 

Total mass of wing, Ib-sec^/ft 

45 

56.85 

26.58 

14.05 

0 

1.458 

2.759 

0.1789 

45 

56.85 

26.58 

14.05 

0 

1.458 

2.759 

0.1644 

Natiiral frequencies, cps: 

Fourth mode 

21 

58 

81 

115' 

21 

52.5 

77 

107 
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TABLE II 


RESULTS OF FLUTTER ANALYSIS AND CORRELATION WITH EXPERIMENT 

(a) Correlation of theory and experiment for model A 
(M = 0.85; P = 0.000787 slugs /cu ft) 



Modes used 
in analysis 

V^, 

fps 

ff, 

cps 


Results of calculations 


1, 2 
1, 3 
1, 2, 3 
1, 2, 3, ^ 

2,1*05 

91*3 

882 

876.5 

73-5 

1*1.5 

1*0.0 

39.8 

0.280 

.1*03 

.1*15 

.1*16 

T*' •V"-r\ ^ 'V' l TYi "n 4" .... 



92I* 

37.9 

0.575 





(b) Results of converged flutter calculations at other Mach numbers and densities 

(Four modes used in analysis) 


M 

P. 

slugs/cu ft 

v^, 

fps 

ff , 
cps 

‘^0,f 

M 

P, 

slugs/cu ft 

fps 

ff, 

cps 

k_ _ 

0,f 

0.85 

0.000504 

.001267 

.0021 

.00526 

1,091* 

720 

618 

554 

! 

59.5 

40.2 

42.5 
45.9 

0.529 

.511 

.630 

.726 

0 

.4 

.7 

.95 


916 

950 

920 

928 

53.5 

52.1 

46.0 

38.0 

0.555 

.515 

.458 

.575 

0 

.7 

.95 

0.000504 

1,129 

I1I55 

1,155 

1,075 

53.0 

51.2 

45.5 

36.9 

0.450 

.414 

.567 

.515 

0 

.4 

*7 

.95 

0.001267 

755 

741 

752 

792 

54.0 

52.6 

46.8 

42.2 

0.673 

.650 

.570 

.488 


(c) Effects of structural modifications 
(M = 0.85; p = 0.000787 slugs/cu ft) 


Case 

Model 

«ij 

(i i J) 

Results of converged 
solutions 

Vf, 

fps 

ff, 

cps 

*^0,f 

1 


0 

797 

38.2 

0.439 

2 


0 

857 

40.0 

.428 

5 ^ 

■ 

Nonzero 

872 

39.0 

.410 


^The model combined the generalized masses of model B and 
the natural frequencies of model A. 
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TABLE III 

GENERALIZED MASSES AND CORRESPONDING CAIOULATED NATURAL 
FREQUENCIES COMPUTED FROM EQUATION (21) FOR MODEL A 



M. . , slugs, for 
ij 

values of j 

of - 


i 

1 

2 

3 

4 

cps 

1 

0.010376 

-0.0015128 

0.0034925 

-0.026760 

20.95 

2 

-.0015128 

.0159856 

-.0159947 

- . 0645656 

57-2 

3 

.003^925 

-.01399^7 

1.245043 

-.0471072 

81.6 

4 

-.026760 

- . 0645656 

-.0471072 

2.897119 

125-7 
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TABLE IV 


MEASURED SLOPES OF NATURAL-VIBRATION MODE SHAPES 

sr 

CALCULATED BY MEANS OF EQUATION (A7)^ 


1 

s 

Semi span 
station, 
percent 

r = 1 

(25 percent 
chord) 

r = 2 

(50 percent 
chord) 

r = 5 

(75 percent 
chord) 

r = 4 
(trailing 
edge) 

■ 

1 

20 

0.280 

0.727 

1.096 

1.350 


2 

35 

1.048 

1.570 

1.870 

2.550 


3 

50 

2.298 

2.964 

3.410 

4.080 

1 

4 

65 

4.001 

4.495 

4.927 

5.284 


5 

80 

5.842 

6.373 

6.755 

6.794 


6 

90 

7-707 

7.957 

7.916 

7.621 


1 

20 

-0.658 

-0.786 

-1.245 

-1.873 


2 

35 

-1.150 

-1.208 

-1.526 

-1.777 


3 

50 

-•993 

-.605 

-.175 

.481 

2 

4 

65 

.992 

1.888 

2.868 

4.012 


5 

80 

6.029 

7.245 

8.685 

10.046 


6 

90 

11.938 

12.558 

12.978 

12.915 


1 

20 

0.891 

0.659 

-0.968 

-5.355 


2 

35 

1.037 

.065 

-1.262 

-5.117 


3 

50 

.783 

-.205 

-1.055 

-1.207 

3 

4 

65 

-.1195 

-.259 

-.116 

1.217 


5 

80 

-.622 

.125 

1.559 

3.656 


6 

90 

-.4505 

1.012 

2.854 

4.987 


1 

20 

0.436 

0.455 

0.667 

1.538 


2 

35 

.303 

.239 

.245 

-.556 


3 

50 

-.591 

-.987 

-1.492 

-5.022 

4 

4 

65 

-1.830 

-1.720 

-1.877 

-2.200 


5 

80 

.373 

1.486 

3.017 

4.788 


6 

90 

7.963 

9.029 

9.196 

8.772 


8L 

Each slope may he obtained by multiplying the corresponding 
value in the table by 0.003225. 
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TABLE V 


MEASURED SLOPES IN PITCH 


a 


(i) 

sr 


OF NATURAL-VIBRATION-MODE SHAPES^ 


i 

s 

Semi span 
station, 
percent 

r = 1 

(25 percent 
chord) 

r = 2 

(50 percent 
chord) 

r = 3 

(75 percent 
chord) 

r = 4 

(trailing 

edge) 


1 

20 

0.1063 

0.204 

0.0968 

-0.166 


2 

35 

.572 

.540 

.484 

.0445 


3 

50 

1.19 

1.16 

.854 

.610 

1 

4 

65 

2.015 

1.85 

1.655 

i. 4 l 6 


5 

80 

2.735 

2.69 

2.62 

2.47 


6 

90 

5.61 

3-355 

3.54 

3.315 


1 

20 

-0.205 

- 0.221 

-0.174 

- 0.164 


2 

35 

-.446 

-.2895 

-.473 

-.755 


3 

50 

.1208 

-.0856 

-.348 

- 1.13 

2 

4 

65 

1.02 

.931 

.499 

-.2455 


5 

80 

5.54 

5.01 

2.98 

2.52 


6 

90 

6.18 

5.80 

5.59 

5.49 


1 

20 

-0.074 

-1.29 

-1.295 

-1.78 


2 

35 

-.556 

-1.713 

-3.05 

- 4.14 

X 

3 

50 

-2.015 

-3.495 

-5.05 

-6.50 

J 

4 

65 

- 4.19 

-5.26 

-6.65 

-7.69 


5 

80 

-6.26 

-6.94 

-7.70 

-8.15 


6 

90 

- 7-55 

-7.58 

- 7-85 

-7.98 


1 

20 

-0.172 

0.067 

0.335 

0.855 


2 

35 

.028 

.084 

.483 

1.426 


3 

50 

-.474 

-.501 

.082 

.479 

4 

4 

65 

.791 

-.483 

-.743 

-1.270 


5 

80 

.790 

.767 

.371 

-.461 


6 

J 

90 

4.515 

4.172 

3.640 

2.937 


a 

Each slope may be obtained by multiplying the corresponding 
value in the table by 0.003225. 
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(a) Control points. 



Flg\ire 1.- Lattices used in evaluation of elements of flutter determi- 
nant for delta semispan wing. 
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Figure 2.- Construction details of delta-wing model. 
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(b) Second mode. Frequency, 52.5 cps. 
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NACA - Langley Field, Va 


Figure Schematic diagram of apparatus used to measure natural -vibration-mode shape 
directions are considered positive as shown; z is considered positive downward. 




